First-principles calculation of the plasmon resonance and of the reflectance spectrum 

of Silver in the GW approximation 
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We show that the position and width of the plasmon resonance in Silver are correctly predicted by 
ab-initio calculations including self-energy effects within the GW approximation. Unlike in simple 
metals and semiconductors, quasiparticle corrections play a key role and are essential to obtain 
Electron Energy Loss in quantitative agreement with the experimental data. The sharp reflectance 
minimum at 3.92 eV, that cannot be reproduced within DFT-LDA, is also well described within 
GW. The present results solve two unsettled drawbacks of linear response calculations for Silver. 

PACS numbers: 71.15.-m, 71.45.Gm, 79.20.-m 
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I. INTRODUCTION 



Although the coexistence of quasiparticles and collec- 
tive excitations in an interacting system is well known, 
our knowledge of their properties and mutual interac- 
tions in real materials is far from being complete. One 
of the most successful approaches to the calculation of 
the quasiparticle (QP) band structure for a wide range 
of materials is the GW methodtl. Starting from a non- 
interacting representation of the system, electrons and 
holes are screened by the surrounding electronic clouds 
created through the Coulomb interaction. This descrip- 
tion has been shown to be rather general in describing 
direct and inverse photoemission spectra, also for compl 
cated systems as reconstructed surfaces and clustersBI 
Noble metals, on the other hand, have been studied only 
very recently. Full quasiparticle calculations have been 
carried out only for Cu obtaining an excellent agree- 
ment with the experimental band structureu. Collective 
modes, i.e. plasmons, occur at energies for which the 
real part of the dielectric function vanishes with a cor- 
responding small imaginary part; they can be observed 
experimentally as sharp peaks in Electron Energy Loss 
Spectra (EELS). A well established technique to calcu- 
late EEL spectra uses the single particle or random phase 
approximation (RPA) for the polarization function, ob- 
tained in terms of Ab-initio energy bands. The latter 
can be calculated within the Kohn Sham (KS)u formu- 
lation of Density Functional TheoryH in the Local Den- 
sity Approximation] (DFT-LDAA or in the Generalized 
Gradient Approximation (GGA)u, when the exchange- 
correlation potential V xc is considered as an approxima- 
tion to the self-energy operator S. This approach com- 
pletely neglects self-energy corrections, that is the differ- 
ence E — V xc , and excitonic effects, whose inclusion has 
been faced only recently. In the case of silicon, Olevano 
and ReiningtJ showed that using the quasiparticle energy 
bands without including excitonic effects the shape of the 
plasmon peak worsens with respect to the experiment. 
Inclusion of self-energy corrections and excitonic effects 



yields a spectrum very similar to DFT-LDA one and to 
the experiment. Ku and EguiluzEJ obtained a correct 
positive dispersion of the plasmon width in K using the 
single particle approximation, with no many-body cor- 
rections beyond DFT-LDA. In these cases many-body ef- 
fects (quasiparticle corrections and/or excitonic effects) 
are not required to describe correctly the experimental 
data. These results agree with the general feeling that ex- 
citonic effects partially cancel self-energy corrections. A 
similar result has been found for CutiEj, where the RPA 
response function calculated without many-body correc- 
tions yields good agreement with the experimental EEL 
and optical spectra. 



In this framework the case of Silver is rather surpris- 
ing: the experimental EELS is dominated by a sharp 
plasmon peak at 3.83 eVc3, whose positipa and width 
are badly reproduced in DFT-LDA RPAEj. In partic- 
ular, a width of about 0.5 eV is obtained within this ap- 
proach, to be contrasted with a much narrower experi- 
mental width (~ 100 mcV). A similar discrepancy occurs 
in the reflectance spectrum, where a very narrow dip at 
3.92 eV is hardly reproduced by DFT-LDA calculations. 
Some papers have recently appearecOE3 which correct 
DFT-LDA results by empirical scissors-operators shifts 
(or similar), meant to better account for the band struc- 
ture. Improved dielectric functions are obtained in this 
way, but no solution of the puzzles mentioned above has 
been given. 



In Sec. II we show that non trivial quasiparticle cor- 
rections on the highly localized d bands strongly modify 
the absorption spectrum threshold. As a consequence, in 
Sec. Ill we show that the plasmon position and width are 
renormalized by quasiparticle effects leading them close 
to the experimental values. Also the reflectance dip at 
3.92 eV turns out to be well described in the same scheme. 



Finally in Sec. IV we summarize our conclusions. 
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II. QUASIPARTICLE EFFECTS ON THE 
ABSORPTION SPECTRUM 



We proceed calculating the band structure of Silver 
within DFT-LDA. The diagonalization of the KS hamil- 
tonian is performed using norm-conserving pseudopoten- 
tials (PPs) and a-^lane waves basial3. The use of soft 
(Martins- TroullierO) PPs allows us to work at full con- 
vergence with a reasonable kinetic energy cutoff (50 Ry.). 
Many-body corrections are added on top of the DFT 
LDA band structure following-thc implementation of the 
GW method described in Refill. The resulting quasipar- 
ticle band structure of Ag at high symmetry points is 
compared with DFT-LDA results in Table @. While the 
deeper energy levels remain mostly unchanged, a down- 
ward shift of about 1.3 eV of the top d bands leads to 
a decrease of the bandwidth, and hence to an excellent 
agreement with experiment. The residual discrepancies 
between the position of the GW bands and the exper- 
imental, results (^0.2eV) is larger than in the case of 
coppera (~0.05eV). This is a consequence of the QP 
renormalization of the Fermi level in silver (~0.6eV), 
not present in copper, which yield the large QP correc- 
tions reported in Table | and an additional error. Differ- 
ently from semiconductors, GW corrections do not act as 
a rigid shift of the whole occupied band structure with 
respect to the empty (conduction) part. QP corrections 
are highly non trivial since even their sign turns out to 
be band/k-point dependent, as we have already shown 
in the case of coppera. 
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TABLE I: Theoretical band widths and band energies for silver, at 
high-symmetry points. GW energies are relative to the QP Fermi 
Level. The striking agreement with the experimental results shows 
that the silver band— structure is very well described at rHie GW 
level. The values in the last column are taken from Ref.Ej where 
spin-orbit splittings have been removed by making degeneracy- 
weighted averages. 



To calculate the EEL spectra, we use the QP band 
structure to evaluate the inverse dielectric function 
e- 1 (w)H 



uj (lu + irj) 



(1) 



where (w) is the interband contribution and lod — 
9.48 eV is the Drude plasma frequency, both calculated 
Ab-initio following the procedure described in Refs.fc3t3. 
The interband RPA dielectric function is given by 



€ ib (Uj) = 1 

n'k- 



E 



A v f d3k 
— 47T nm / 77 

i^°Jbz (2tt) 3 

I 2 

■ q|e~ iq " r |nk) 



In' 



k-q 



q 



(2) 



where < /„ (k) < 2 represents the occupation number 
summed over spin components. 77 = 0.1 eV has been used 
in the present calculations. The q — > limit of Eq. (^) 
has been done including the effects, of the pseudopoten- 
tial non-locality, as done in Ref.113. The GW optical- 
transition energies (-E n ,k — ^«',k) are obtained fitting the 
QP corrections calculated at 29 k-points of a regular 
grid in the irreducible Brillouin zone (BZ). Separate fit- 
ting curves (shown in Fig. |l|) have been used for selected 
band pairs (n, n') and k regions in the BZ, in order to 
correctly reproduce the energy dependence of the quasi- 
particle corrections. 
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FIG. 1: GW optical-transition energies (circles) plotted as a func- 
tion of DFT-LDA ones for a regular grid of 29 k— points in the 
irreducible Brillouin zone. The dashed line corresponds to vanish- 
ing GW corrections while full lines are our fitted curves. They are 
used to extend the GW corrections to the optical transition ener- 
gies of a larger k-points grid. GW corrections shift upward the 
electron-hole energies, except for some transitions between 4 and 
5eV (as discussed in the text). In the inset GW corrections are 
shown versus DFT-LDA energies. 



From Fig. it is evident that GW corrections are 
highly non-trivial (ranging from —0.75 to 0.5 eV for 
filled states, and from to 0.75 eV for empty states, 
as shown in the inset). GW optical-transition en- 
ergies are generally above the dashed-line, represent- 
ing the condition of vanishing quasiparticle corrections 



E n > tk ) = fj5° jk --Ev )k ), and lead to upward 
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shifts of the electron-hole energies. However between 4 
and 5 eV some electron-hole energy shifts are negative. 
These correspond to transitions close to the L point from 
the Fermi level to the 7th (empty) band. 

In Fig. U the GW interband contribution e ib (u>) is com- 
pared with the DFT-LDA result. In the upper panel all 
the transitions except those involving the k-points near 
the L point are included: GW corrections shift the whole 
spectra to higher energies, with the threshold energy oc- 
curring at ~ 4eV, in agreement with experiment. In 
panel (b) we consider the transitions (sketched in the in- 
set) not included in panel (a) . At difference with panel 
(a), GW shifts the DFT-LDA spectrum toward lower en- 
ergies. In conclusion, weak interband transitions start at 
the threshold in panel (b) (3.36 eV), while the main on- 
set of e ib {ui) is that of panel (a), ~ 4eV. The resulting 
(total) dielectric function is shown in the inset of Fig. ^. 



through the relation 




FIG. 2: Interband absorption spectrum e ib (ui) of Ag calculated 
using the quasiparticle GW band structure. Theoretical spectra do 
not contain the intraband contribution. In both panels and insets: 
solid line, GW; dashed liae, DFT-LDA. Only in panel (a): boxes, 
full experimental e (oj)Ex In panel (a) all the transitions except 
those involving the k-points near the L point are included. These 
are considered separately in panel (b). In the insets the DFT- 
LDA band structure is compared with the GW result, and the 
arrows indicate the most important transitions involved near the 
absorption threshold. In panel (b) quasiparticle corrections shift 
some transition energies at the L point below the main threshold 
of panel (a). These are responsible for the plasmon damping. 



III. THE RENORMALIZED PLASMON 
FREQUENCY 



The two different effects of GW corrections on e ib (u) 
shown in Fig. || are crucial in determining the properties 
of the plasmon resonance whose frequency, ui p , is defined 



0. 



(3) 



lu p is, in general, complex since decay into electron- 
hole pairs gives to the plasmon a width proportional to 
e ib (uj p ). Now, the sharp onset of e ib (oj) at ~ 4eV (see 



inset of Fig. |3j) is responsible for a solution of Eq. (|3j) at 
3.56 eV, just below the main interband threshold. Within 
DFT-LDA, instead, the interband threshold of e ib (u>) at 
~ 3 eV remains below u> p , resulting in a strongly damped 
plasmon peak as shown in Fig. [?] (dashed line). From 
Fig. ||, panel (b) , we see that the only transitions con- 
tributing to the plasmon width, proportional to e ib (u) p ), 
are those from the Fermi surface to the first conduc- 
tion band at the L point. These decay channels are in 
agreement with many experimental resultsEj that, using 
temperature and alloying techniques, have singled out 
the electron-hole transitions responsible for the plasmon 
damping. 
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FIG. 3: Electron Energy Loss Spectrum (EELS) of Sihiea. Solid 
line: GW. Dashed line: DFT-LDA. Boxes: experiment^. The 
non trivial quasiparticle GW corrections improve considerably the 
DFT-LDA plasmon peak, yielding a striking agreement with the 
experiment. The GW dielectric function is shown in the inset. 



GW-corrected EELS is compared with DFT-LDA and 
experimental results in Fig. |: the plasmon peak, un- 
derestimated in intensity and position in DFT-LDA, is 
shifted toward higher energies and strongly enhanced by 
GW corrections, in striking agreement with experiment. 
In Fig. ^ an overestimation of the absorption spectrum 
with respect to ±he experiment is observed (similar to 
that found in CuE3) . This is reflected in an underestima- 
tion of the loss intensity close to 8 eV, and might be due 
to the neglection of QP renormalizajtion factors in Eq. (||) 
and to excitonic dynamical effectscd. Despite this partial 
agreement in the 8eV region, the present Ab-initio GW 
calculation is able to reproduce correctly the plasmon 
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resonance. This resonance can be interpreted as a col- 
lective (Drude like) motion of electrons in the partially 
filled band. However, according to Eq. |], lo p does not co- 
incide with the bare Drude frequency cJd, the difference 
arising from the screening of the electron-electron inter- 
action by virtual interband transitions. The plasmon res- 
onance, although blueshifted with respect to DFT-LDA, 
remains below the main interband threshold, but overlaps 
the weak low-energy tail of interband transitions (see 
Fig. ||, Frame (b)) acquiring a small -yet finite- width. 
The highly non-trivial QP shifts of interband transitions 
are crucial to obtain this result. 

The polarization of the "medium" where the plasmon 
oscillates (the d electrons) is hence important to deter- 
mine its energy and width. This polarization is absent 
in the homogeneous electron gas because there are no 
localized d orbitals and no interband transitions; it is 
weak in semiconductors (like Si) , because interband tran- 
sitions occur at energies far from that of the plasma res- 
onance. The same polarization effect is present, but de- 
structive in copper due to the lower onset of interband 
transitions. EELS peaks occur above this onset and are 
therefore strongly broadened. In conclusion, the delicate 
interplay of plasmon-frequency renormalization with the 
shift of the interband-transition onset, both due to QP 
corrections, may yield (in Silver) or may not yield (in 
Copper) a sharp plasmon resonance. 

Another important quantity is the reflectance, R (lu) = 
(\N (w)-l\/\N (uj) + 1|) 2 , where N is the complex re- 
fraction index defined by [N (oj)] 2 = e(uj). In Fig. |] 
we compare the GW R (col with the DFT-LDA one, and 
with experimental resultsE3. The latter shows a very nar- 
row dip at 3.92 eV, close to the plasmon frequency, arising 
from the zero-reflectance point <x>o, defined as e (wo) = !• 
Again, the width and depth of this reflectance dip are 
related to the imaginary part of e(cu). GW corrections 
make c^o to occur below the main onset of interband 
transitions, and hence produce a very narrow and deep 
reflectance minimum. Here the agreement between GW 
results and experiments for the intensity and width of the 
dip at 3.92 eV is even more striking than in the EELS. 

This result is of great importance for the optical and 
EELS properties of Ag surfaces. Very recent calculations 
of Reflectance Anisotropy Spectra (RAS) within DFT- 
LDA for Ag(110)E3 were not able to reproduce quantita- 
tively a sharp dip observed experimentally. This feature, 
at an energy uj r = 3.8 eV, has been assigned to a bulk 
resonance, arising in the RAS spectrum when e (u> r ) = 1. 



Hence, it is the same occurring in the reflectance spec- 
trum of Fig. f| Its width and shape are strongly related 
to the reflectance dip of Fig. |^, and therefore they need 
GW corrections to be well reproduced. 
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FIG. 4: Reflectivity spectrum of Silvej. Solid line: GW. Dashed 
line: DFT— LDA. Boxes: experimentEj. The experimental sharp 
dip at 3.92 eV is correctly reproduced by GW, with a substantial 
improvement on the DFT-LDA spectrum. 



IV. CONCLUSIONS 

We have performed a calculation of the EEL and reflec- 
tivity spectra of silver within the RPA approach, using 
the quasiparticle band structure calculated within the 
Ab-initio GW method. We have shown that the pecu- 
liar, well known plasmon peak in the EELS and the deep 
reflectivity minimum observed experimentally, are quan- 
titatively described by theory, for the first time without 
requiring adjustable parameters. Our theoretical calcu- 
lations do not contradict previous results for semiconduc- 
tors and simple metals (where many-body effects are less 
important) if the role played by d orbitals is correctly 
interpreted. 
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